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摘 要 : 植被 生物 量 反 映 了 生态 系统 获取 能 量 的 能 力 ,分 析 其 分 布 特征 对 了 解 生态 系统 结构 和 功能 具有 十 分 重要 


的 意义 。 传 统 的 反 演 植被 地 上 生物 量 的 方法 往往 由 于 样本 的 缺少 ， 
影像 灰 度 值 和 植被 因子 等 10 个 因子 ) 进行 训练 ， 


块 样 地 进行 验证 ,结果 表明 :ELM 反 演 植被 地 上 生物 量 , 可 以 获得 较 高 的 精度 ,模型 预测 结果 与 实测 结果 的 曲线 拟 


高 。 本 文选 用 ELM 对 105 块 实测 样本 的 遥感 因子 (T 


头 及 影响 因子 的 不 确定 性 而 导致 预 佑 精度 不 
用 余下 34 


合 决定 系数 尽 达 0.89。 此 外 ,对 2010 一 2015 年 玛 纳 斯 河流 域 的 植被 地 上 生物 量 进行 反 演 ,认为 流域 内 上 游 山 区 
生物 量 大 部 分 较为 稳定 ,中游 平 原 区 生物 量 呈 现 增加 趋势 ,下 游 荒 漠 区 生物 量 则 呈现 退化 趋势 。 
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生物 量 (biomass ) 泛 指 单位 面积 上 所 有 生物 有 
BLPRU- 38 CU 。 研 究 不 同 尺度 下 植物 地 上 生物 量 
的 空间 分 布 格局 ,不仅 对 了 解 植物 物种 多 样 性 .生产 
力 和 群落 动态 具有 指示 作用 ,而 且 对 了 解 植 物种 类 
分 布 与 外 界 环境 的 关系 也 具有 重要 的 参考 价值 ” 。 
通常 线性 模型 并 不 能 很 好 地 表达 生物 量 和 光学 反射 
率 之 间 复 杂 的 非 线 性 关系 中 ,而 光学 影像 的 NDVI 
主要 反映 植被 冠 层 信 息 , 不 能 全 面 地 获取 植被 的 立 
体 信 息 中 。 此 外 ,由 于 地 表 参 数 的 多 变性 和 复杂 
性 ,传统 以 经 验 、 半 经 验 模型 为 主 的 生物 量 反 演 方法 
往往 只 对 生物 量 和 后 向 散射 系数 的 关系 作 各 种 线性 
或 非 线性 的 拟 合 5 E ,忽略 了 各 种 统计 函数 无 法 表 
达 所 有 的 不 确定 性 。 而 人 工 神经 网 络 则 非常 适用 于 
表达 这 种 非 线 性 的 映射 关系 ,可 以 充分 逼近 任意 复 
杂 的 非 线 性 关系 ,学 习 和 自 适 应 不 知道 或 不 确定 的 
系统 ,有 很 强 的 鲁 棒 性 和 容错 性 ”" 。 所 以 ,用 人 工 
神经 网 络 模型 估算 植被 生物 量 可 以 说 是 开辟 了 一 条 
方便 快捷 可 靠 的 途径 和 方法 。 

考虑 到 温带 干旱 地 区 不 同 土地 利用 类 型 的 生物 
量 信息 仍然 缺乏 ,本 文选 取 干 旱 区 典型 流域 一 一 玛 
纳 斯 河流 域 作为 研究 区 ,根据 TM 影像 光谱 信息 以 
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及 NDVI 等 植被 指数 ,结合 地 面 调查 数据 对 研究 区 
植被 地 上 生物 量 建立 神经 网 络 反 演 模型 。 通 过 分 析 
玛 纳 斯 河流 域 2010 年 与 2015 年 植被 地 上 生物 量 的 
空间 分 布 规律 ,为 合理 利用 干旱 区 资源 ,保护 其 生态 
平衡 提供 重要 理论 依据 。 


1 材料 与 方法 


1.1 研究 区 概况 

玛 纳 斯 河流 域 地 处 天 山北 荔 中 段 , 准 嘲 尔 例 地 
南 缘 , 南 北 纵 路 纬度 2°05' (85? 01' ~ 86*32' E, 
43?27' -45?21' N) 。 北 部 与 古 尔 班 通 古 特 沙漠 相 
接 , 南 部 为 山地 ,中 游 的 绿洲 区 由 冲积 洪 积 鹿 .泉水 
溢出 带 、 冲 积 平原 三角洲、 湖滨 平原 构成 5 2] 。 玛 
纳 斯 河流 域 行政 区 域 包括 石河 子 垦 区 , 沙 湾 县 . 玛 纳 
斯 县 .新 湖 总 场 ( 图 1) ,流域 内 地 势 从 南 向 北 逐 渐 降 
低 , 地 貌 呈 典型 的 山 倪 系统 结构 ,依次 为 山地 一 山 前 
倾斜 平原 一 沙漠 。 由 于 地 貌 、 气 候 、 水 文 条 件 . 植 被 
具有 明显 的 垂直 地 带 性 分 布 " -中 ,土地 利用 类 型 复 
杂 多 样 (以 草地 、 林 地 、 耕 地、 荒漠 为 主 ) ,生物 群落 存 
在 明显 的 区 域 性 和 差异 性 ,使 得 本 研究 具有 典型 性 。 
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A1 玛 纳 斯 河流 域 分 布 示 意图 


Fig.1 Distribution of the Manas River Basin 


1.2 遥感 数据 来 源 

本 研究 所 需 遥 感 数据 有 :GD 基础 地 理 数 据 , 包 
括 流域 边界 .各 县 市 边界 、 河 流 等 ,来 自 于 新 疆 测绘 
局 ;@) 土地 利用 数据 ,包括 2010 年 和 2015 年 的 土 
地 利用 矢量 数据 ,该 数据 来 源 于 30 m 分 辨 率 土地 覆 
盖 数 据 集 ,土地 利用 类 型 按 林 地 、 草 地、 耕地 ,湿地 、 
人 工 建设 用 地 ,其 他 (主要 包括 苔 葡 / 地 衣 、 裸 岩 、 戈 
BE RE 沙漠 .盐碱地 .冰川 /永久 积 雪 等 ,以 沙漠 为 
主 ) 等 6 大 类 进行 划分 ,成 图 比例 尺 1: 250 000; 3) 
TM 遥感 影像 ,来源 于 地 理 空 间 数据 云 30 m 分 辩 率 
的 Landsat 系列 数据 (条 带 号 143/144, 行 编号 29/ 
30) ;由 DEM 数字 高 程 数据 ,来 源 于 地 理 空间 数据 
云 GCDEMDEM 30 m 分 辩 率 数字 高 程 数据 。 以 上 数 
据 均 经 过 新 疆 生 产 建 设 兵 团 环境 监测 中 心 站 验证 ， 
符合 要 求 ( 注 : 所 有 数据 在 ArcMap10.2.2 及 ENVI 
5.3 支持 下 ,通过 格式 转换 .边界 裁剪 .几何 校正 、 空 
间 又 加 等 处 理 , 进 行 数据 统计 分 析 , 所 有 数据 坐标 系 
统一 为 Albers Conical Equal Area 投影 ,基准 面 为 
WGS -1984 ,中央 经 线 为 105^E) 
1.3 不 同 植被 样 地 地 上 生物 量 数据 来 源 
1.3.1 森林 生态 系统 地 上 生物 量 地 面 观测 

(1) 样 地 设置 :林地 固定 样 地 面积 为 10 000 
m , 样 方 2 个 ,面积 为 30 m x30 m。 林 下 植被 样 方 3 
个 ,面积 为 5 m x 5 m, 具 体 以 对 角 线 中 心 点 及 两 边 
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各 距 15 m 处 为 中 心 点 ,布设 3 个 5 mx5 m 的 区 域 。 
参照 4 林地 分 类 》(LYZT1812 - 2009 ) 和 森林 资源 规 
划 设 计 调 查 主要 技术 规定 (国家 林业 局 2003 年 ) 进 
行 观测 。 

(2) 观测 内 容 : 主 要 包括 林木 种 类 株数、 胸径、 
郁 闭 度 、 叶 面积 指数 .生物 量 等 。 

(3) 生物 量 计算 :森林 生态 系统 地 上 生物 量 观 
测 分 为 立木 地 上 生物 量 和 和 冠 层 下 部 活体 植被 地 上 生 
物 量 观测 ,立木 地 上 生物 量 与 冠 层 下 部 活体 植被 地 
上 生物 量 之 和 即 为 样 方 生物 量 。 立 木 地 上 生物 量 观 
W :通过 样 方 内 所 有 林木 进行 测量 ,获取 其 树 高 . 胸 
径 等 地 面 观 测 数据 ,依据 相对 生长 方程 计算 ,对 所 有 
立木 生物 量 求 取 平 均值 并 除 以 样 方面 积 ,获取 1 m? 
面积 的 立木 生物 量 。 冠 层 下 部 活体 植被 地 上 生物 量 
观测 :在 样 方 内 ,随机 选择 3 个 5 m x5 m 的 区 域 ,分 
别 收集 其 中 全 部 地 上 植被 , 称 量 鲜 重 ,并 从 中 抽取 不 
少 于 5% 的 样品 ,105 C 下 烘 干 称 干 重 ,获取 植株 含 
水 量 ,进而 获得 实测 的 地 上 生物 量 ,计算 3 个 区 域 平 
均值 并 除 以 样 方面 积 ,作为 冠 层 下 部 1 m 面积 的 生 
物 量 。 

1.3.2 草地 生态 系统 地 上 生物 量 地 面 观 测 

(1) 样 地 设置 :草地 固定 样 地 面积 为 10 000 
m , 取 7 一 9 月 植物 生长 期 观测 。 观 测 采 用 现场 调 
查 法 ,参照 《青海 省 草地 资源 调查 技术 规范 》 
( DB63/ F209 - 1994 ) 和 《草地 旱 鼠 预测 预报 技术 规 
程 )(DB63/T331 - 1999) 设 置 样 地 。 

(2) 观测 内 容 :主要 包括 草地 羡 度 、 叶 面积 指 
数 、 生 物 量 等 。 

(3) 生物 量 计算 :按照 草地 生态 系统 样 方 布设 
规则 布设 9 个 1m x1 m 样 方 ,将 样 方 内 植物 地 面 以 
上 所 有 绿色 部 分 取 下 并 做 标记 。 称 量 鲜 重 后 ,65 C 
烘 干 称 量 干 重 并 记录 。 计 算 9 个 样 方 活 体 生物 量 之 
和 ,将 9 个 值 取 平 均 ,得 到 1 m 面积 的 生物 量 。 
1.3.3 农田 生态 系统 地 上 生物 量 地 面 观 测 

(1) 样 地 设置 ;农田 固定 样 地 面积 为 10 000 
m^ ,农田 样 方面 积 为 1 mxl m, 参照 《 基 本 农田 划 
定 技术 规程 》(TDZT1032 -2011 ) 进行 观测 。 

(2) 观测 内 容 : 主 要 包括 作物 群体 株 高 . 叶 面 积 
指数 .生物 量 等 。 

(3) 生物 量 计算 :采用 样 方 法 , 样 方 大 小 为 
1 m x1 m, 至 少 需 9 个 重复 。 农 田 样 地 选择 要 远离 
树木 、 田 间或 建筑 物 , 样 方 选择 距离 路 边 田 十 、 沟 边 
等 3 m 以 上 。 在 作物 成 熟 后 收获 前 的 晴天 实施 野外 
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调查 ,采集 作物 地 上 部 分 全 植株 体 。 按 照 农田 生态 
系统 样 方 布设 规则 布设 9 个 1 m x1 m 样 方 ,将 样 方 
内 所 有 地 上 植株 剪 取 并 标记 。 返 回 实 验 室 后 ,将 作 
物 植 株 分 为 叶 、 茎 .籽粒 3 部 分 分 别称 量 鲜 重 并 记 
录 , 在 80 烘箱 中 烘 干 至 恒 重 ,记录 干 重 数据 。 计 
算 9 个 样 方 烘 干 生物 量 并 记录 ,将 9 个 值 取 平均 ,得 
到 1 m 面积 的 生物 量 。 
1.4 数据 处 理 
1.4.1 逐 感 数据 及 处 理 ” 建 模 所 采用 的 遥感 数据 
JJ Landsat TM 影像 ,时 相 为 2012 年 8 月 。 实 地 测量 
的 植被 生物 量 数据 时 间 为 2012 年 8 月 中 旬 ,影像 时 
相 较 一 致 ,植被 变化 较 小 ,影像 云 量 为 0.4% ,能 够 
清楚 地 识别 植被 因子 和 地 物 。 

首先 对 Landsat TM 影像 在 ENVI 5.3 平台 下 进 
行 几何 校正 .数据 投影 变化 、 辐 射 定 标 `, 大气 校正 以 
及 研究 区 的 裁剪 ,然后 在 ENVI classic 平台 下 导入 
野外 测 得 的 样 地 点 坐标 ,并 提取 影像 7 波段 的 像 元 
值 。 之 后 在 AveGIS 平台 下 将 野外 测 得 的 样 地 点 坐 
标 导入 ,并 与 波段 运算 后 的 影像 进行 琶 加 ,得 到 样 地 
点 的 比值 植被 指数 (RVI) 、 归 一 化 植被 指数 ( ND- 
VI) 、 差 值 植被 指数 (DVI) 、 绿 度 (GREEN ) 7 m E 
(FG) 和 SQRT(RVI1) 等 ,最 后 利用 ArcGIS 的 3D An- 
alyst 工具 对 DEM 数据 的 坡度 (slope) 、 坡 向 (aspect) 
因子 进行 提取 。 
1.4.2 极限 学 习 机 由 于 传统 的 神经 网 络 训练 时 
间 过 长 .计算 开销 大 、 对 初始 参数 设置 敏感 \ 易 陷入 
局 部 极 值 、 收 敛 速 度 慢 其 至 不 能 收 钱 、 过 拟 合 问 
题 “。 因 此 ,本 文 用 单 隐藏 层 前 馈 神经 网 络 一 一 极 
限 学 习 机 ( extreme learning machine, ELM ) jr iUe 
对 研究 区 植被 地 上 生物 量 建立 神经 网 络 模型 ,并进 
行 误差 分 析 。 

图 2 为 一 个 单 隐 藏 层 前 馈 神 经 网 络 基本 原理 
图 ,网 络 中 及 个 输入 节点 ,L 个 隐 层 节点 ,WM 个 输 
出 节点 ;w; 是 第 i 个 隐 层 节点 的 输入 权重 问 量 ;b; 是 


图 2 极限 学 习 机 基本 原理 示意 图 


Fig.2 Basic principle diagram of the extreme learning machine 


第 i 个 隐 层 节点 的 偏差 (影响 因子 );pB; 为 第 i 个 隐 
层 节 点 输出 到 第 j 个 输出 节点 的 权重 值 ;g (x) 是 激 
活 函数 。 

对 于 任意 给 定 的 输入 数据 集 (X,Y) | (x%i ,yi)， 
(x2 ,72) (33,y3) (Xm,ynm)| FEP x; = [a ,Xi ,Xi 
ix]: 是 第 i 个 样本 的 输入 ,y; = [ya ,yo Yos 
Yin 为 数据 集中 第 i 个 样本 的 输出 。 第 i 个 样本 输 
入 到 极限 学 习 机 (ELM) 中 ,第 j 个 隐藏 层 节 点 的 输 
入 如 式 (1): 


WiX; = (oj 二 OP tst +wp) Xa 


二 OOjiXil tOpXp5 tOgXa 十 cb OX, (1) 


经 过 第 7 个 隐 层 节点 处 理 后 ,第 7 个 隐 层 节点 输 
出 g(oyx; * bi) o L 个 隐 层 节点 的 输出 矩阵 为 : 


h(x;) = ig Cox; *tb,) gCosx; *b,) 
g(waxi +b) ge( wx, * b,)| (2) 
根据 隐藏 层 输出 矩阵 五 和 输出 权重 矩阵 BB 可 
以 通 近 输出 矩阵 了 : 
Bi 


p- 
Y=h(x:) B= h(x) B, 


Ba 
= |g(wix; +b) g(omx; +b) gCosx; +b) 
Pu WEN Bin 
g Cox, t b,)] : 
Pu WE Bin 
2(0,,0,,0,:--0,) (3) 


极限 学 习 机 (ELM ) 的 激活 函数 无 线 可 微 , 可 以 
无 限 允 近 输出 值 , 即 零 误差 的 估计 输出 值 。 那 么 存 


在 合适 的 obb EEE, = 0。 所 以 , 利 
用 最 小 二 乘法 , 求 得 输出 矩阵 B: 

B-H'Y (4) 
APH” SES HUBER RE ZH. AI UAE E ; 


M Y NKIRE. HIRERE B 时 ,极限 学 习 机 的 
学 习 过 程 就 结束 了 。 
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1.4.3 ”基于 神经 网 络 的 生物 量 反 演 模型 第 一 步 : 
相关 性 分 析 。 对 提取 的 因子 (TMI ,TM2 ,TM3 ,TM4、 
TMS TM7 RVI NDVI, DVI, GREEN, FG, SQRT 
(RVI) „aspect ,slope) 进行 相关 性 分 析 , 最 终 选 取 了 
TM1 ,TM2,TM3, TMA , TM5,TM7 和 NDVI、RVI、DVI 
以 及 SQRTCRVI) 等 10 个 相关 性 较 高 的 因子 作为 神 
经 网 络 输入 层 自 变量 ,实测 生物 量 为 输出 层 变量 。 
第 二 步 : 归 一 化 处 理 。 为 了 将 不 同 量 纲 .不 同 数 
量 级 的 原始 数据 放 在 一 起 比较 ,通常 需要 对 原始 数 
据 以 及 代表 统计 指标 进行 变换 处 理 。 本 文选 用 规格 
化 变换 法 "…” , 即 从 原始 数据 的 矩阵 中 找到 每 个 变量 
的 最 大 值 和 最 小 值 ,然后 用 原始 数据 和 最 小 值 的 差 
除 以 最 大 值 和 最 小 值 的 差 , 即 按 下 面 的 公式 对 每 个 
数据 进行 转换 : 
X, - (X, -a)/(b -a) (5) 


式 中 :X 为 原始 数据 ;X; 为 归 一 化 处 理 后 的 数据 ;oa 为 
最 小 值 ;b 为 最 大 值 ; 其 中 X; 的 值 都 介 于 [0 ,1 之 间 。 

第 三 步 : 结 构 选 择 与 参数 设置 。 将 实测 植被 生 
物 量 数据 (包括 森林 、 农 田 与 草地 ) 分 为 两 部 分 ,一 
部 分 进行 训练 , 男 一 部 分 进行 验证 ,通过 调试 隐藏 神 
经 元 数目 , 找 出 最 优 参数 。 


2 结果 与 分 析 


2.1 基于 神经 网 络 模型 的 生物 量 反 演 结果 

随机 选取 105 块 样 地 的 相关 遥感 因子 作为 神经 
网 络 训练 模型 的 输入 矩阵 ,与 之 对 应 的 实测 生物 量 
样 地 数据 作为 神经 网 络 训练 模型 的 参照 值 ,进行 神 
经 网 络 模型 的 训练 。 并 用 余下 的 34 块 样 地 相关 遥 
感 因子 作为 极限 学 习 机 的 预测 模型 的 输入 矩阵 ,用 
来 估 测 实际 生物 量 。 通 过 实际 生物 量 与 模型 预 佑 生 
物 量 进行 比 对 ,验证 神经 网 络 模型 的 精准 度 。 图 3 


隐藏 神经 元 数 =20 


隐藏 神经 元 数 =30 
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为 不 同 隐藏 神经 元 数 的 预测 生物 量 与 实测 生物 量 的 
对 应 关系 。 

从 图 3 可 以 看 出 :模型 预 佑 的 生物 量 与 实测 生 
物 量 相关 性 较 高 , 当 隐 藏 的 神经 元 数 较 小 或 较 大 时 , 
神经 网 络 模型 的 精度 都 会 有 所 降低 。 为 了 清楚 直观 
地 表达 神经 网 络 模拟 生物 量 的 精度 ,笔者 将 神经 网 
络 预测 值 与 实测 值 进行 数据 拟 合 。 拟 合 函 数 有 : 指 
ZI PR ( Custom Equation ) , [8i FE HH pj ( Fourier) 、 高 
斯 函数 (Gauss) 插值 函数 (Interpolant) .多项式 函数 
( Polynomial ) , X£ pj Zt ( Power ) I IE 5% PR 2C ( Sum of 
Sin) 。 通 过 比较 ,选择 拟 合 效果 最 好 的 函数 值 作为 该 
隐藏 神经 元 下 神经 网 络 预测 模型 的 拟 合 函数 ,与 之 对 
应 的 相关 系数 作为 预测 模型 精度 的 评判 标准 。 表 1 
为 不 同 隐藏 神经 元 数 的 预测 值 与 实测 值 的 拟 合 结果 。 

如 表 1 所 示 : 在 模型 的 预测 值 与 实际 值 拟 合 中 ， 
(EENT pL ži ( Fourier) , 3 ir PAZ ( Gauss ) 使 用 的 较 
多 , 且 拟 合 效果 较 好 ,其 他 函数 如 指数 函数 ( Custom 
Equation) ) 插值 函数 (Interpolant) , Z Ji 5X pK% ( Pol- 
ynomial ) FAE PA ( Power ) 拟 合 效果 较 差 ,不 能 较 好 
地 模拟 出 预测 值 与 实际 值 的 关系 ;隐藏 神经 元 数 = 
50 时 ,模型 的 预测 值 与 实际 值 拟 合 效果 最 好 ,曲线 
拟 合 决定 系数 R i 0. 89, 且 预测 值 与 实测 值 的 线 
性 相关 系数 最 高 , 达 0. 67; 当 隐 藏 神 经 元 数 =20 及 
60 时 ,曲线 拟 合 的 决定 系数 R 均 为 0.79。 由 此 可 
充分 说 明 用 神经 网 络 模 型 反 演 生 物 量 是 可 行 的 。 

通过 对 比 , 选 择 隐藏 神经 元 数 = 50 作为 反 演 遥 
感 生物 量 的 模型 参数 。 图 4 为 隐藏 神经 元 数 =50 
时 模型 模拟 的 34 块 样 地 的 预测 值 与 实测 值 的 拟 合 
效果 。 

从 图 4 可 以 看 出 : 当 隐藏 神经 元 数 50 ERE, BS 
型 反 演 的 耕地 生物 量 与 草地 生物 量 较为 精准 ,但 模 
拟 林地 生物 量 时 ,部 分 样 地 模拟 值 要 低 于 实测 值 。 


隐藏 神经 元 数 =40 隐藏 神经 元 数 =50 


一 一 隐藏 神经 元 数 =60 一 一 隐藏 神经 元 数 =70 一 一 一 实测 值 


10 000 
8 000 
6 000 
4 000 
2 000 


生物 量 /(g - m?) 


样品 号 
注 :林地 样本 号 1 ~2,10 ~12,15 ~21,25 ~26; 草 地 样本 号 3 ~6,9,29 ~30; 耕 地 样本 号 7 ~8,13 ~14,22 ~24,27 -28,31 ~34。 下 同 。 
图 3 神经 网 络 模型 反 演 生 物 量 结 


Fig.3 Results of inversion biomass of the neural network model 
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表 1 不 同 隐藏 神经 元 数 的 预测 值 与 实测 值 的 拟 合 结果 
Tab.1 Fitting results of the predicted and measured values of different hidden neurons 
R-square Adjusted R-square RMSE Corrcoef 拟 合 函 数 
隐藏 神经 元 数 =20 0.792 9 0.755 9 1 179 0.495 6 4l; FR nF PR ( Fourier) 
隐藏 神经 元 数 =30 0.1619 0.078 1 2 292 0.233 1 在 里 叶 函 数 ( Fourier) 
隐藏 神经 元 数 =40 0.422 5 0.133 7 2222 0.308 8 高 斯 函数 (Gauss) 
隐藏 神经 元 数 =50 0. 890 8 0.8103 1 040 0.667 6 正弦 函数 (Sum of Sin) 
隐藏 神 经 元 数 =60 0.7917 0.687 6 1 334 0.391 8 高 斯 函数 ( Gauss ) 
隐藏 神经 元 数 =70 0.2510 0.1173 2 243 0.186 3 在 里 叶 函 数 ( Fourier) 
注 :R-square 代表 曲线 拟 合 决定 系数 ,Adjusted R-squared 代表 改良 后 的 曲线 拟 合 决定 系数 ,RMSE 代表 标准 差 , Corrcoef 代表 预测 值 与 实测 值 的 


线性 相关 系数 。 


生物 量 /(g * m?) 


图 4 隐藏 神经 元 数 =50 时 4 


O 隐藏 神经 元 数 =50 


EI Scil ft 


011121314 151617 18 19 20 2122232425 2627 28 29 30 31 323334 
样本 号 
E 物 量 实测 值 与 预测 值 拟 合 结果 


Fig.4 Fitting results of the predicted and measured values when the number of hidden neurons was 50 


分 析 其 原因 发 现 :用 神经 网 络 佑 测 遥 感 影像 生物 量 ， 
虽然 拟 合 效果 较 好 ,但 由 于 遥感 影像 的 影响 因子 较 
复杂 , 且 样 地 上 植株 分 布 无 规律 性 ,导致 反 演 结果 有 
偏差 ,因此 ,用 神经 网 络 反 演 植 被 生物 量 还 存在 一 定 


缺陷 。 
2.2 玛 纳 斯 河流 域 不 同 土地 利用 类 型 的 生物 量 分 
布 特征 


为 了 探究 不 同 土地 利用 类 型 生物 量 的 空间 分 布 
特征 ,在 ArcMap10. 2.2 平台 下 对 2010 年 与 2015 年 
30 m 分 辨 率 的 China Cover 土地 利用 数据 集 进行 林 
地 草地、 耕地 提取 ,将 顶 格 转 面 ,分 别 加 载 神经 网 络 
模型 反 演 的 2010 年 与 2015 年 8 月 植被 地 上 生物 量 
数据 ,最 终 得 到 2010 年 与 2015 年 林地 .草地 和 耕地 
三 种 土地 类 型 的 植被 地 上 生物 量 分 布 图 (图 5)。 

图 5 所 示 :2010 年 玛 纳 斯 河流 域 各 土地 利用 类 


地 总 体 平均 生物 量 为 2 634.47 g m^, 

至 2015 年 , 玛 纳 斯 河流 域 各 土地 利用 类 型 中 ， 
草地 仍 占 主导 地 位 , 占 地 面积 为 12 841.47 km? , 平 
均 生 物 量 为 239. 01 g - m? , H6 2010 年 无 明显 变 
化 ;耕地 面积 为 7 953.18 km ,面积 较 2010 年 增加 
T 504.99 km? ,平均 生物 量 为 1 341.90 g-m”, Æ 
物 量变 化 较 明 显 , 可 能 是 种 植 结构 发 生 了 变化 ;2015 
年 林地 面积 为 1 572. 90 km ,总 平均 生物 量 为 
4470.21 g * m 一, 虽然 平均 生物 量 相 较 2010 年 有 
明显 增加 ,但 从 林地 生物 量 空间 分 布 来 看 ,在 玛 纳 斯 
河流 域 的 下 游荡 漠 区 有 和 较 明 显 的 林地 消失 ,退化 面 
HIK 657. 26 km ,人 类 应 该 加 强 对 林地 的 保护 , 合 
理 开 发 利用 林地 。 

总 的 来 说 ,2010 一 2015 年 玛 纳 斯 河流 域 上 游 林地 
和 草地 生物 量 有 小 幅度 增加 趋势 ,中游 耕地 呈现 生物 


型 中 , 占 地 面积 最 大 的 为 草地 ,主要 分 布 在 玛 纳 斯 河 
流域 上 游 山区 ,面积 达 12 737. 16 km ,平均 生物 量 
为 281.96 g * m”, 占 研究 区 总 面积 的 37.41% ;耕地 
面积 为 7 448. 192 km ,主要 分 布 在 玛 纳 斯 河流 域 中 
游 平原 区 ,平均 生物 量 为 919. 68 g. m ^ Matri EU 
2 230.16 km ,主要 分 布 在 研究 区 的 上 游 山 区 及 下 游 
沙漠 区 ,其 中 上 游 山 区 平均 生物 量 为 8 746.53 g - 
m“ ,下游 沙漠 区 平均 生物 量 仅 为 414.37 g - m^, 


量 增加 趋势 ,而 下 游 沙 漠 区 生物 量 呈现 退化 趋势 。 
结论 与 展望 


本 文 以 TM 影像 像 元 灰 度 值 和 植被 因子 作为 模 
型 输入 层 变量 ,利用 单 隐 藏 层 前 馈 神 经 网 络 一 一 极 
限 学 习 机 (ELM ) 方 法 对 玛 纳 斯 河流 域 地 上 生物 量 
进行 反 演 。 结 果 表 明 : 虽 然 ELM 是 通过 “黑箱 ” 操 
作 , 难 以 归纳 出 指导 性 规律 ,但 该 算法 能 够 克服 传统 
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(c) 2010 年 耕地 


| 生物 量 /(g * m?) 
EN 0 ~ 500 


生物 量 /(g * m?) 
m0 ~ 200 


E 200 ~ 400 E 500 ~ 1 000 

= 400 ~ 600 = 1 000 ~ 1 500 
Em 600 ~ 800 E 1 500 ~ 2 000 
E800 ~ 1 000 = 2 000 ~ 10 000 


(f) 2015 年 耕地 N 
0 1020 40 60km 


生物 量 /(g * m?) 


EB 0-200 EH 0-500 

ES 200 ~ 400 = 500 ~ 1 000 
= 400 ~ 600 æ 1 000 ~ 1 500 
ES 600 ~ 800 ES 1 500 ~ 2 000 
= 800 ~ 1000 = 2 000 ~ 10 000 


注 : 由 于 林地 草地、 耕地 所 含 植被 地 上 生物 量 不 同 , 故 分 类 标准 也 不 同 。 


5 玛 纳 斯 河流 域 生物 量 分 布 


Fig.5 Distribution of biomass in the Manas River Basin 


梯度 算法 (如 BP 算 法 ) 局 部 最 小 化 以 及 过 度 拟 合 的 
问题 ,最 大 限度 地 利用 遥感 影像 样本 集 的 先 验 知识 ， 
在 反 演 植 被 地 上 生物 量 时 可 以 获得 较 高 的 精度 。 
通过 对 2010 一 2015 年 玛 纳 斯 河流 域 植 被 地 上 
生物 量 的 反 演 ,结果 表明 : 玛 纳 斯 河流 域 草地 生物 量 
一 般 为 0~200 gm 一, 主要 分 布 在 上 游 山 区 和 下 
游 沙 漠 区 , 近 5 a 来 无 明显 变化 ;林地 生物 量 跨度 较 
大 ,上 游 山 区 林地 生物 量 平均 值 高 达 8 746.53 g - 
m ^ ,而 下 游 沙 漠 区 林地 生物 量 平均 值 仅 为 414. 37 
g:m , 且 近 5 a 来 沙漠 区 退化 较 明 显 , 林 地 退化 面 
HIK 657.26 km? ;中 游 平原 区 多 为 耕地 ,2010 一 2015 
年 生物 量 平均 值 变 化 较 显 著 ,5 a 间 生 物 量 平 均值 增 
加 了 395.22 g. m ,可 能 是 种 植 结构 变化 所 致 。 
极限 学 习 机 (ELM ) 学 习 速 度 快 且 不 必 过 分 注 
意 数据 集 的 大 小 ,但 在 反 演 实测 生物 量 时 会 遇 到 两 
个 常见 的 问题 。 第 一 ,激活 函数 的 选择 问题 ,由 于 激 
活 函数 的 选择 往往 是 从 竺 处理 问 题 中 提取 先 验 的 知 
识 结合 进 极限 学 习 机 (ELM) 中 ,而 实际 应 用 中 获取 
先 验 知 识 往往 比较 困难 。 第 二 ,变量 间 相 关 性 问题 ， 


即 求解 广义 逆 和 矩阵 时 无 法 避免 输入 相关 数据 造成 的 
病态 解 或 权 值 过 大 等 问题 。 这 可 能 是 神经 网 络 模型 
在 预测 生物 量 时 存在 缺陷 的 原因 。 
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Abstract: Aboveground biomass reflects the capability of ecosystems to obtain energy. Analysis on the spatial dis- 
tribution pattern is of great significance for understanding the structure and function of ecosystems. The accuracy of 
inverting aboveground biomass with the conventional approach is low due to the lack of samples and the uncertainty 
of impact factors. In this study , Extreme Learning Machine ( ELM) was used to train the remote sensing factors of 
105 samples which included seven-band pixel values of TM image and vegetation factors , and the remaining 34 sam- 
ples were used for verification. The results confirmed that ELM approach could invert vegetable aboveground bio- 
mass with a higher accuracy ,and its determination coefficient of curve fitting reached 0. 89. In addition ,the inver- 
sion of vegetation aboveground biomass in the Manas River Basin from 2010 to 2015 found that the biomass was rel- 
atively stable in the upper area of the Manas River Basin , was in an increase trend in the middle plains and was in 
a deterioration trend in the downstream desert. 
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